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Abstract 

We argue that the approximate "Casimir scaling" of the string tensions of higher- 
representation Wilson loops is an effect due to the finite thickness of center vortex config- 
urations. It is shown, in the context of a simple model of the Z 2 vortex core, how vortex 
condensation in Yang-Mills theory can account for both Casimir scaling in intermediate 
size loops, and color-screening in larger loops. An implication of our model is that the 
deviations from exact Casimir scaling, which tend to grow with loop size, become much 
more pronounced as the dimensionality of the group representation increases. 



1 Introduction 



There is increasing numerical evidence |]-f|] supporting the Center Vortex Theory of quark 
confinement [§-0, which was put forward in the late 1970's. Briefly, a center vortex is 
a topological field configuration which is line-like (in D=3 dimensions) or surface-like (in 
D=4 dimensions) having some finite thickness. Creation of a center vortex can be regarded, 
outside the line-like or surface-like "core," as a discontinuous gauge transformation of the 
background, with a discontinuity associated with the gauge group center. Creation of a 
center vortex linked to a Wilson loop, in the fundamental representation of SU(N), has 
the effect of multiplying the Wilson loop by an element of the gauge group center, i.e. 

W(C) -> e t27Tn/N W(C) n = 1,2,..., AT- 1 (1) 

The vortex theory, in essence, states that the area law for Wilson loops is due to quantum 
fluctuations in the number of center vortices linking the loop. 

Paradoxically, this emphasis on the center of the gauge group can be viewed both as 
a vital strength of the theory, and also as a fatal weakness. Both aspects are apparent 
when we consider the force between static quarks in an SU(N) gauge theory, whose color 
charge lies in the adjoint representation. The QCD vacuum will not tolerate a linear 
potential between adjoint quarks over an infinite range; this is simply because adjoint 
color charges can be screened by gluons. Asymptotically, the force between adjoint quarks 
must drop to zero, and this is exactly what happens in the center vortex theory. The 
adjoint representation transforms trivially under the group center; adjoint Wilson loops are 
unaffected by center vortices, unless the core of the vortex happens to overlap the perimeter 
of the loop. As a result, large loops have only perimeter falloff, and the force between 
adjoint quarks vanishes asymptotically. The argument extends to any color representation 
which transforms trivially under the center of the gauge group. The fact that center 
vortices make such a clear distinction between those color charges which should be confined, 
and those charges which should not, is one of the most attractive features of the theory. 

The fatal weakness aspect was first pointed out in ref. [|K|. Consider the large- N limit, 
which has the factorization property < AB >=< A >< B >, where A and B are any two 
gauge-invariant operators. Then 

< W A (C) >=< W F (C)W F {C) >=< W F (C) > 2 (2) 

where Wa,f{C) denotes, respectively, Wilson loops in the adjoint and fundamental repre- 
sentations. An immediate consequence is that confinement of fundamental representation 
quarks implies confinement in the adjoint representation, with string tension a a = 2o-p. 
This is possible because color screening by gluons is a 1/A^ 2 suppressed process, so at large 
N the vacuum can support an adjoint string. But adjoint loops are insensitive to center 
vortices, as noted above. The apparent conclusion is that center vortices cannot be the 
confinement mechanism at large N. 

Even more troubling is the fact that the existence of an adjoint string tension is not 
peculiar to large- N. Many numerical experiments in 577(2) and SU(3) lattice gauge theory 
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have shown that flux-tubes form, and a linear potential is established, between quarks in the 
adjoint (and higher) representations [fTXf| . The string tension is representation-dependent, 



and appears to be roughly proportional to the eigenvalue of the quadratic Casimir operator 
of the representation. Thus, for an SU(2) gauge theory 

4 

Oj ~ -]{] + 1M/2 (3) 

where <jj is the string tension in representation j. The region where this relation is valid, 
from the onset of confinement to the onset of color-screening, we call the "Casimir-scaling 
regime" Hl2| , |l3|j . Of course, the color charge of higher- representation quarks is eventually 
screened by gluons, and the force between quarks then depends only on the transformation 
properties of the representation with respect to the gauge group center; i.e. on the "n-ality" 
of the representation. Asymptotically, for an SU{2) gauge group, 

a . = [ °i/2 J = half-integer ,^ 
3 1 j = integer 

Color screening, although it must occur for adjoint quarks at sufficiently large separation, 
is very difficult to observe in numerical simulations.^] Existing Monte Carlo studies of the 
QCD string have mainly probed the Casimir scaling regime. 

In short, there is a linear potential between adjoint quarks in some finite range of 
distances, with approximate Casimir scaling of the string tensions, even at iV = 2. Casimir 
scaling should become exact, in a region extending from the confinement scale to infinity, 
in the N —>■ oo limit. Yet according to the center vortex theory, it would appear that string 
formation between adjoint quarks is impossible, at any N, at any distance scale. This has 
always seemed to us a good reason for discarding the vortex theory. 

But suppose — and in our view the numerical evidence ]I|-0 is becoming persuasive 
— that center vortices really are the confining configurations, at least for quarks in the 
fundamental representation. Then either there is some other mechanism for inducing a 
linear potential between adjoint quarks, or else there must be a loophole in the "fatal 
weakness" arguments. The first alternative does not seem very economic, and in any case 
we have no insight, at present, in that direction. We will concentrate instead on the second 
possibility because there is, in fact, one possible loophole. 



2 Thick Vortices and the Loop Perimeter 

The statement that adjoint loops are unaffected by center vortices contains one slight 
caveat: They are unaffected unless the vortex core somewhere overlaps the perimeter of 
the loop. At first sight this caveat seems irrelevant; for large loops the effect can only 
contribute to the perimeter falloff. But suppose that the vortex thickness is actually 

1 At least, this is difficult at zero temperature. Color screening has been observed in certain finite- 
temperature studies, cf. Miiller et al. in ref. jLl[ . 
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quite large: on the order of, and perhaps exceeding, the typical diameters of low-lying 
hadrons. What would be the effect of vortices on Wilson loops whose area is smaller than, 
or comparable to, the vortex cross section? We will study this question in the context of 
a simple model of vortex/perimeter overlaps, mainly in SU(2) lattice gauge theory. 

In D=3 dimensions, or on a constant-time hypersurface in D=4 dimensions, a vortex is 
a closed tube of magnetic flux. For simplicity we consider planar and, in D=4 dimensions, 
spacelike Wilson loops. If the vortex is linked to a Wilson loop, with winding number = 1, 
then the vortex pierces the minimal area of the loop an odd number of times. If the vortex 
is not linked to the loop, it either does not pierce the minimal area at all, or pierces it an 
even number of times. If, for the moment, we ignore the finite radii of the vortex tubes, 
then the effect of vortices on a Wilson loop is simple: For every instance where the minimal 
surface is pierced by a center vortex, insert a center element —I somewhere along the loop, 
i.e. 

W{C) = Ti[UU....U] -> Tr[UU..{-I)..U] (5) 

In principle we should place the —I at the point of discontinuity of the gauge transformation 
which creates the vortex. However, since —I commutes with everything, the placement 
is arbitrary; and this is related to the fact that the Dirac sheet of a center vortex can be 
moved about by gauge transformations. 

Denote by / the probability that any given plaquette on the lattice is "pierced" by a 
vortex; i.e. a line running through the center of the vortex tube intersects the plaquette. 
The area law for Wilson loops is then trivially derived from the assumption that these 
probabilities, for plaquettes in a plane, are uncorrelated.0 In that case, one has 

<w(c)> = n{(i-/)+/(-i)}<^ (c)> 

= exp[-a{C)A] < W {C) > (6) 

where the string tension is 

<7 = -ln(l-2/) (7) 

and where < Wq(C) > is the expectation value of the loop with the constraint that no 
vortices pierce the minimal area. The quantity < Wq{C) > can be (and has been) computed 
from lattice Monte Carlo, cf. ref. In those computations, it is found that W (C) does 

not have an area law falloff. 

By the same argument, the string tension for loops in any j = half-integer representation 
is the same as for j = 1/2, while the string tension for j = integer vanishes. Of course the 
argument is too simple in a number of respects, e.g. there is likely to be some short-range 
correlation between the / probabilities of nearby plaquettes. This point, however, is not 
crucial to the discussion. What is more important is that we have ignored the finite radii 
of the vortices. Equation (||) is only true if the core of the vortex, where it crosses the 
plane of the loop, is entirely contained in the minimal area of the loop. If the core overlaps 
the perimeter of the loop, then eq. (|) cannot be quite right. 

2 Some related ideas are found in ref. 0]. 
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What is needed is a full-fledged theory of center vortices, perhaps something along 
the lines of the old Copenhagen vacuum 0, which would explain how eq. @ should be 
modified when the vortex core is not entirely enclosed within the loop. In lieu of that, we 
will just consider a simplified picture in which the center element —I in (|5[) is replaced by 
a group element G, which interpolates smoothly from —I, if the core is contained entirely 
with the loop, to +/, if the core is entirely exterior. Our assumption, for Wilson loops in 
any group representation j, is the following: 

Assumption 1 

The effect of creating a center vortex piercing the minimal area of a Wilson loop may 
be represented by the insertion of a unitary matrix G at some point along the loop 

W{C) = Tt[UU....U] -> Tt[UU..G..U] (8) 

where 

— * 

G(x, S) = exp[iac(x)n ■ L] 

= S exp[ia c (x)L3]S'^ (9) 

The Li are group generators in representation j, n is a unit 3- vector, and S is an 
SU (2) group element in the j-representation. 

The parameter otc{x) £ [0, 2n] depends on what fraction of the vortex core is enclosed 
by the loop; thus it depends on both the shape of the loop C, and the position x of the 
center of the vortex core, relative to the perimeter, in the plane of the loop. If the core is 
entirely enclosed by the loop, then = 2tt, conversely, if the core is entirely outside 

the minimal area of the loop, then ac{x) = 0. In an abelian theory, this first assumption 
would be completely correct, where Tr[G] would be the value of the loop for a vortex 
created on a classical vacuum background. In a non-abelian theory the assumption might 
be quantitatively correct for expectation values (i.e. averaging over group orientations S 
and over small quantum fluctuations around the vortex background); this would be 
quite sufficient for our purposes. 

Generalizing a little further, if we create some number m of vortices in the loop, centered 
at positions xi, x%, x m , then 

W{C) - W[C;{xi,Si}] 

= Tr[U..UG(x a ,S a )U...UG(x b ,S b )U...UG(x p ,S p ) U) (10) 

where a,b, ...,p is some permutation of 12..., m. We now make the second assumption of 
our model: 

Assumption 2 

The probabilities / that plaquettes in the minimal area are pierced by vortices are un- 
correlated. The random group orientations associated with Si are also uncorrelated, 
and should be averaged. 
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These two assumptions define our model. They are, no doubt, an oversimplification of 
the effects of vortex thickness, but we believe they at least provide a plausible picture of 
those effects. 

According to the second assumption, we are justified in averaging independently every 
G(x a , S a ) over orientations in the group manifold specified by S a . This is easily seen to 
give 

G(a) = J dS S expliaL^S^ = Qj [a] hj+i 
1 1 j 

= 27TT Trexp ^ aL ^ = 2tTT ^ cos ( am ) 

J J 771= — ) 



1 „ r . „ , 1_ ' 

1 

m=—j 

(11) 



sin[(2j + l)f 



(2j + 1) sin[f ] 

where Ik is the k x k unit matrix. Then W[C, {xi, Si}}, averaged over all Si, becomes 

W[C; {x t }] = {l[g j [ac(x i )}}Tr[UU.....U} (12) 

i 

Next, take the expectation value of Tr[UU.....U] for configurations U with the constraint 
that no vortices pierce the loop, denoted < Wq(C) >. Eq. fljjP then goes to 

< WiC-ix,}} >= ([[QiMxi)]} < W Q (C) > (13) 

i 

The last step is to sum over the number and position of vortices piercing the plane of the 
loop C, weighted by the appropriate probability factors, and we find 

<W(C)> = Y[{(l-f) + fg j [ac(x)}}<W (C)> 



exp 



£M(i-/) + /$M*)]} 



< W (C) > 



= exp[-a(C)A] < W (C) > (14) 

where 

®c = ~EM(1 - /) + m<*c{x)]} (15) 

A X 

The product and sum over positions x run over all plaquettes in the plane of the loop. The 
reason for not restricting x to lie strictly within the minimal area of the loop is, again, 
because the vortex core is finite. Denote the radius of the vortex core by R c . If the center 
of the core lies outside the loop, but at a distance less than R c from the perimeter, then 
it can still overlap the perimeter. One can restrict the sum to run over x G A', where A' 
includes all plaquettes inside the minimal area in the plane of the loop, as well as plaquettes 
in the plane outside the perimeter, but within a distance R c of the loop. 

Now <7c is n °t exactly a string tension, because it depends, via ac(x), on the shape (and 
the area) of the loop. If, however, there is some region where oc changes only slowly with 
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area, then the potential will rise approximately linearly. In particular, consider the limit of 
very large loops. In that case, almost every vortex which affects the loop is entirely enclosed 
by the loop, and for these vortices otc(x) ~ 2ir. Only those vortices near the perimeter have 
ctc(x) different from 2tt, and as the loop becomes very large this is a negligible fraction of 
the total; in particular, A' /A « 1. This means that oc is an area- independent constant 
for large loops, and it can be seen from eq. (p"T|) and ([14]) that 

f — ln(l — 2f) j = half-integer 

°c = { n • , 16 

[ (J j = integer 

which is the correct representation-dependence of the asymptotic string tension. 

Next, let us consider the case where / << 1, which is certainly true in the lattice theory 
at weak coupling, and also small or medium size loops, where otc{x) is also typically small. 
For small a, we have from (|TT|) 

Si[a]«l- yj(j + l) (17) 
Then, making an expansion of the logarithm in (|TUj) and applying (|TTD, 

A 



°c = /~£(i-$M*)]) 

xeA' 



or just 

where o? c is an average value 



vc = ^ c3 ( 3 + l) (19) 



-2 1 

'V — t 



nr- ^£«c(z) (20) 

We see that ac, for small / and small loops, is proportional to the eigenvalue of the 
quadratic Casimir operator. 

This result can be readily extended to any SU(N) group. In the general case there are 
N — 1 types of center vortices, corresponding to the N — 1 phase factors of eq. flU). To the 
n-th type, we associate probability /„ to pierce a plaquette, and a group factor 

G[x, S] = Sexp[ia%{x) ■ H]S^ (21) 

where the {Hi, i = 1,..,N — 1} are the generators spanning the Cartan sub-algebra. 3 
Following the same steps as before 

<w(c)> = n(i- 

x { n=\ ) 

3 It is possible that only vortices with the smallest magnitude of center flux have substantial probability; 
i.e. fi = /jv-1 is finite, all other /„ are negligible. This is a dynamical issue which we cannot resolve here, 
and so we consider the general case that includes all possible /„. 
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Q r [a\ = -j-Trexp[ia • #] (22) 

with d r the dimension of representation r. Vortices of type n and type N — n have phase 
factors in (|1|) which are complex conjugates of one another; they may be regarded as the 
same type of vortex but with magnetic flux pointing in opposite directions, so that 

/„ = /*_„ and g T [a n c {x)}=g;[a%- n {x)} (23) 

and therefore 

<W{C)> = n{l- E/n(l-Re^[a^)])| 

X I n=l J 

a c = -^£ln|l- £ / n (l-R*£ P [^(z)])J (24) 

Expanding the logarithm to leading order in f n , expanding G r [®\ to leading order in a, and 
using the identity 

1 



one finds 

°c = T 1 2^ i-nr(-r) • o, (.r) 




C( 2) (26) 



where C^ 2 - 1 is the eigenvalue of the quadratic Casimir operator of the SU(N) group in 
representation r. 

The results ( p"9D and (p^) might be termed "Casimir proportionality," since the non- 
perturbative part of the interquark potential, which is due to vortices, is proportional to 
the quadratic Casimir of SU(N) for small loops. But this does not yet imply Casimir 
scaling of string tensions. The parameters Sq(x) depend on loop size, and there is no 
particular reason to suppose that oc is constant in the adjoint representation or, equiv- 
alently, that the adjoint potential is linear in some range. Even if the adjoint potential 
were approximately linear in some interval, it is not obvious that the string tension for 
the fundamental representation, in the same range of distances, would have reached its 
asymptotic value. To study this issue, we will return to the SU(2) example. 



3 Linear Potentials and ctc{x) 

It may be possible to measure «c(i) in computer simulations, by the methods introduced 
in refs. In the meantime, it is worthwhile to ask whether there exists some reasonable 

ansatz for ac(x) which would lead to both Casimir proportionality and linear potentials 
in some region. 

To set things up, let us consider a long rectangular RxT loop in the x — t plane, with 
T » R, in group representation j. The time-extension T is huge but fixed, so we can 
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characterize loops C just by the width R. Let x denote the x-coordinate of the center of 
a vortex, where it pierces the x — t plane. From the previous discussion, the inter-quark 
potential induced by vortices is easily seen to be 



Yi(R) = - E M(i -/) + /$[«*(*»)]} (27) 

n=— oo 

where x n = n + | (the choice of x n comes from the fact that the vortex centers lie in the 
dual lattice, piercing the middle of plaquettes). The problem is to find some reasonable 
ansatz for or(x). 

Suppose the time-like sides of the loop are at x = and x = R. Then, to guide the 
search for an ansatz, there are a few simple conditions that 0£r(x) must satisfy: 

1. Vortices which pierce the plane far outside the loop don't affect the loop. This means 

that for fixed R, as x — > ±oo, we must have (Xr(x) — > 0. 

2. If the vortex core is entirely contained within the loop, then olr(x) = 2ir. This 
translates as follows: Let x be inside the loop, and d be the distance from x to the 
nearest of the time-like sides. Then it must also be the case that 0£r(x) — > 2tt as 
d — > oo. 

3. As R — > 0, the percentage of any vortex core which is contained inside the loop must 

also go to zero. Thus an(x) — > as R — > 0. 

There are an infinite number of functional forms which would meet these conditions, but 
a simple 2-parameter (a, b) ansatz is the following: First define 

f \ [ x — R for \R — x\ < \x\ , nnS 

U\ x ) — \ c r> (28) 

yv ' 1 —x for \R — x\ > \x\ v ' 

whose magnitude is the distance of the vortex center to the nearest timelike side of the 
loop, taken negative if the vortex center is inside the loop, and positive outside. Then 
choose 

1 — tanh ay(x) + — 



n /,'l.r) = 7T l — tann ^ay{x) + —J (29) 

which fulfills all three requirements. 

Figure 1 shows the potentials for the j = ~, 1, § representations, for the choice of 
parameters / = 0.1, a = 0.05, b = 4., in the range R G [1, 12] Note that the fundamental 
and adjoint potentials are roughly linear in a range from 5 or 6 to 12 lattice spacings. 
Figure 2 plots the ratios Vx(R) /Vi/ 2 (R) , an d V 3 / 2 (R) /Vi/2(R). As expected, these ratios 
start out at the ratios of the corresponding Casimirs, 8/3 and 5 respectively. It can be seen 
that in this interval the adjoint/fundamental ratio drops only slowly, from 8/3 to about 2, 

4 Strictly speaking, R takes on only integer values in the lattice formulation, but we have plotted Vj(R) 
over the continuous interval. 



9 



cc 
> 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 




Vortex-Induced Potentials: j=1/2, 1,3/2 



'fund' 
'adjoint' 
'3|2' 




Figure 1: Interquark potential V(i2) induced by center vortices, according to the model 
discussed in the text, for quark charges in the in the j — |, 1, § representations. 
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Figure 2: Ratios of V-i(R) (upper curve) and V\(R) (lower curve) to the potential Vi(R) 
of fundamental representation quark charges. 
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Vortex-Induced Potentials: j=1/2, 1, 3/2 
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Figure 3: Same as Figure 1, for the range R G [1, 100]. 



while in the same interval the (j = 3/2) /fundamental ratio drops more precipitously, from 
5 to about 2.5. Figure 3 is again a plot of all three potentials, but this time over the range 
[1, 100]. There are two things to notice in this last figure. First, we see color screening set in 
as expected, around R — 25 — 30 lattice spacings. The adjoint potential goes to a constant, 
while the j = 3/2 potential closely parallels the fundamental potential. Secondly — and 
this was not expected — it can be seen that the the slope of the fundamental potential, in 
the Casimir scaling region between 6 to 12 lattice spacings, is very close to its asymptotic 
value. 

We have therefore found, in this rather simple model of the vortex core, the kind of 
Casimir scaling which is seen in Monte Carlo simulations. A natural question is to what 
extent this scaling depends on a very special choice of parameters. The answer is that 
Casimir scaling, although it is not found for any choice of parameters, is generic in a large 
region of the parameter space. For example, the Casimir scaling region of Fig. 3 can be 
scaled up by any factor F simply by setting a — > a/F, b — * bF. It is also quite easy to 
think up other functional forms for an(x) which satisy the above three conditions (e.g. 
olr{x) = (3(x) — f3(x — R), with (3{x) — > ±n in the limits x — > ±oo). The possibilities 
are, of course, infinite. Some functions work better than others, but the existence of an 
approximate Casimir scaling region of some finite extent seems to be fairly common. The 
crucial ingredient is the thickness of the vortices, which would be on the order oil/ a for 
the choice of ocr(x) in eq. (p9|). What we find is that the thickness of the vortices must 
be quite large; larger, in fact, than the Casimir region itself, in order to see approximate 
Casimir scaling of the adjoint string tension. 
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4 Conclusions 



We have presented a scenario whereby the Casimir scaling of higher-representation string 
tensions is explained in terms of the finite thickness of center vortices. We do not claim 
to have proven that vortex thickness is the origin of Casimir scaling, but this explanation 
now appears to be very plausible, particularly if center vortices turn out to be the true 
QCD confinement mechanism. 

Numerical tests of our scenario are in order. If center vortices give rise to an adjoint 
string tension, and if an adjoint loop is evaluated only in those configurations where no 
vortex links the loop, then the string tension should vanish. This was found to be the case 
for the fundamental string tension, in refs. 0,0) and should be testable for the adjoint 
string, at least in principle, by the methods explained in those articles. It may also be 
possible to calculate ctn(x) from Monte Carlo simulations of fundamental loops, use that 
information in eq. (0) to derive the adjoint potential, and compare the derived adjoint 
potential with the corresponding Monte Carlo data. 

For every set of parameters in which we see Casimir scaling in our simple model, 
the deviations from exact Casimir scaling are much greater for the j = 3/2 potential as 
compared to those of the j = 1 potential. It would therefore be very interesting to compare 
Fig. 2 with the actual Monte Carlo data. Calculation of the j = 3/2 potential in the scaling 
region of SU(2) gauge theory, by lattice Monte Carlo techniques, is a computationally 
intensive problem, but we believe it is feasible. 

Finally, there is the question raised years ago in ref. [TD|: how can center vortices explain 
confinement at large N, where the Casimir scaling region extends to infinity? We have 
seen, in the scenario outlined here, that the size of the Casimir scaling region depends on 
the thickness of center vortices. Already in SU(2) gauge theory, the existence of a Casimir 
regime implies that the diameter of the vortices substantially exceeds the separation length 
at which the heavy quark potential begins to grow linearly (i.e. the "confinement scale"). 
We therefore expect that as N increases, the diameter of vortices relative to the confinement 
scale will also slowly increase, probably as log(iV). At iV = oo the vortex core would be 
infinite in extent, and the discontinuous gauge transformation exterior to the core would 
be pushed off to infinity. It is not clear that the term "center vortex" remains a useful 
description of the relevant configurations in this limit. 



Note Added - We have recently learned that a related proposal, namely, that center 
vortices might lead to a (breakable) string potential between massive gluons, was advanced 
by Cornwall |ff5| in a conference proceedings from 1983. We thank Mike Cornwall for 



bringing this reference to our attention. 
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